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Thermal noise of optical reference cavities sets a fundamental limit to the frequency instability 
of ultra-stable lasers. Using Levin's formulation of the fluctuation-dissipation theorem we correct 
the analytical estimate for the spacer contribution given by Numata et al. [I]. For detailed analysis 
finite-element calculations of the thermal noise focusing on the spacer geometry, support structure 
and the usage of different materials have been carried out. We find that the increased dissipation 
close to the contact area between spacer and mirrors can contribute significantly to the thermal 
noise. From an estimate of the support structure contribution we give guidelines for a low-noise 
mounting of the cavity. For mixed-material cavities we show that the thermal expansion can be 
compensated without deteriorating the thermal noise. 



I. INTRODUCTION 

Ultrastable optical cavities have become a standard 
tool for stabilizing laser systems needed e.g. for high- 
resolution spectroscopy, optical clocks [3J [3] , optical mi- 
crowave generation [JHS] and coherent optical frequency 
transfer [7HPJ • State-of-the-art cavity-stabilized laser sys- 
tems show linewidths below 1 Hz and fractional frequency 
instabilities between 10 -16 and 10~ 15 at one second jTDJ- 

□2]. 

The fractional length stability 6L/L of ultra-stable 
cavities is limited by inevitable thermal noise in the cav- 
ity materials and a world-wide effort is ongoing to re- 
duce this limitation. An efficient way of calculating ther- 
mal noise has been proposed by Levin |18j , applying the 
fluctuation-dissipation theorem. Following his so-called 
"direct approach" for optical cavities, a set of equations 
was given by Numata et al. |T] for estimating the noise 
contribution of spacer, substrates and coatings. 

The equations give a good estimate of the thermal 
noise arising from the mirrors which have been identified 
as the dominant source of noise for cavities with mirror 
substrates made of the widely used ultra low expansion 
glass (ULE). In state-of-the-art optical cavities however, 
substrate materials of higher mechanical quality are be- 
ing used [P7 1H9H2T] . For these optical cavities a fractional 
instability on a level of a few 10 -16 has been shown. At 
this level the spacer contribution is non-negligible and 
has to be calculated by numerical modeling. In this pub- 
lication we describe the mechanisms of dissipation in an 
optical cavity with focus on the spacer contribution, in- 
cluding the support structure, the usage of different ma- 
terials and the spacer geometry. We discuss current cav- 
ity designs with respect to their thermal noise budget 
and, where possible, will give construction guidelines for 
future cavity designs. 

In section [TT| we introduce the theoretical framework 



used for the simulations. We apply the simulation to a 
typical cavity design in section |III A| and discuss devi- 
ations from the analytical estimate for the spacer con- 
tribution. In section |III B| we show that also the cavity 
support can lead to non-negligible contribution to the 
spacer thermal noise. In the final section [ill C| we focus 
on the thermal noise floor of mixed-material cavities. 



II. THEORETICAL FRAMEWORK 

Thermal noise in laser mirrors has been investigated in 
great detail with respect to gravitational wave detection 
[151 1221- 28 . A variety of thermal noise sources have been 
identified so far |29) . Brownian thermal noise being the 
most dominant one. An effective way to estimate thermal 
noise has been first proposed by Levin [IB] , following the 
fluctuation-dissipation theorem as introduced by Callen 
and Welton [30]. The effect of Brownian motion ther- 
mal noise on the length stability of ultra-stable optical 
resonators has been pointed out by Numata et al. [T] 
following Levin's so called "direct approach" This 
concept is illustrated for clearness in the following. 

For optical resonators the generalized coordinate de- 
scribing the length fluctuations is the averaged distance 
between the two cavity mirrors x — x\ — X2 as probed by 
the laser field. According to the fluctuation-dissipation 
theorem, the spectral density S x (f) of the thermal fluc- 
tuations of x is calculated by applying a corresponding 
conjugate force with amplitude Fq to each mirror surface 
leading to 



SM) 



2k B T W di: 



7T 2 / 2 



F 2 



(1) 



where T denotes the temperature and ks the Boltzmann 
constant. The term Wdi ss denotes the time- averaged dis- 
sipated power in the system when an oscillatory force 
with amplitude F and frequency / is applied. To calcu- 
late the thermal noise for a homogeneously distributed 
internal loss the term Wdi ss is then expressed by 
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where U denotes the maximum elastic strain energy and 
and 4> the loss angle of the system. The total dissi- 
pated power and thus the noise of the cavity length x 
can be described as the sum of the contributions from 
spacer S% (/), two substrates Si (/) and two coatings 
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The power spectral density of length fluctuations S x (f) 
can easily be converted to fractional frequency fluctua- 
tions S y (f) — S x (f)/L 2 . The instability of the fluctua- 
tions in the time-domain is conventionally characterized 
by the Allan-deviation a y [31] of the fractional frequency 
fluctuations y. The 1// noise of equation (|3| (flicker fre- 
quency noise) leads to a constant Allan deviation 
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In a typical reference cavity the differential mirror dis- 
placement is probed by a Gaussian laser beam with a 
1/e 2 beam radius w. Consequently a pressure distribu- 
tion, 



(5) 



has to be applied to the two mirrors with opposite sign 
to drive the excitation. 

If the mirror is treated as an infinite half space the 
contribution of the substrate is given by [T51 H2] 
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where E denotes Young's modulus and a Poisson's ratio 
of the substrate material. Bondu et al. [22] have later 
refined the formula for a finite-size mirror, revealing that 
a long "bar-shaped" mirror substrate with a thickness 
larger than the mirror radius should be preferred to a 
"gong-shaped" substrate. However, for typical cavity ge- 
ometries the beam waist is much smaller than the mirror 
dimensions and therefore the mirror can be approximated 
in good accuracy by a infinite half-space. 

The coating contribution to Brownian motion thermal 
noise is typically estimated by treating the coating as a 
thin layer of thickness d c t on the substrate's front face 
[26), l27l [32] . Assuming a homogeneous loss angle </> ct and 
similar elasticity of coating and substrate as in the case of 
Ti2C>5/Si02 coatings on ULE or fused silica, the coating 
contribution reduces to 
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For materials with large Young's modulus such as silicon 
or sapphire the thermal noise is dominated by the loss 
angle for strains perpendicular to the surface which is 
difficult to access experimentally (see discussion in [35] 
for details). 

As a rough estimate for the spacer contribution to the 
cavity noise the thermal fluctuations of the length of a 
cylinder, averaged over the whole front face area, are 
calculated. A cylinder of length L, radius i? sp and central 
bore radius r sp is assumed. Levin's direct approach then 
demands a pressure p = F /A sp uniformly distributed 
across its front faces with area A sp , leading to an elastic 
energy of 
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According to equations and ^ the thermal fluctua- 
tions of the spacer length L averaged over the full cross- 
section of the front faces are 
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Note, that this estimate differs from the result given by 
Numata et al. in two ways. Firstly, equation ^ takes 
into account that the front face area of the spacer is re- 
duced by the area of the central bore. Secondly, our result 
is a factor of 3/2 larger than the result by Numata et al. 
To derive their equation the authors used the formula for 
the position fluctuations of one end of a free, elastic bar 
[33] . This approach is widely employed for single mirrors 
in interferometers for gravitational wave astronomy. To 
calculate the length fluctuations of the spacer of a Fabry- 
Perot interferometer the authors simply added the fluc- 
tuations from both sides which assumes that those are 
uncorrelated. According to Levin's direct approach the 
calculation of the length fluctuation requires to apply op- 
posite forces simultaneously to both ends. Consequently 
this approach automatically takes into account correla- 
tions from both ends which were neglected in pQ. 

According to the model described by equation ([9]), the 
noise can be minimized by distributing the stress on a 
large cross section. However, for mirrors smaller than the 
diameter of the spacer, this approximation is not valid 
and for a reliable calculation of the spacer contribution 
of the thermal noise S x sp ^ finite element methods (FEM) 
have to be applied as described in the next section. 



III. SIMULATIONS 

For the calculations presented in the following chapter, 
the software package COMSOL [SJ has been used. The 
general cavity geometry is depicted in Figure [T] 

The calculations have been carried out in axial sym- 
metry. Because of mirror symmetry with respect to 
the cavity mid plane, half a cavity was simulated. Fig- 
ure [2] shows the calculated cavity deformation under the 
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FIG. 1: Sketch of the cavity model with dimensions used for 
the FEM simulations in this publication 



pressure distribution of equation ^ assuming the pa- 
rameters given in Table [i] For illustration purposes a 
beam waist of 2 mm has been chosen. We quote our 
results in terms of elastic strain energies. For ULE a 
strain energy of 1 nJ corresponds to a length fluctuation 
Sl = 8.58 • 10~ 35 m 2 /Hz at 1 Hz and a corresponding 
Allan deviation of a y — 1.09 • 10 -16 at 1 second for a 
cavity length of 100 mm. 

The shape of the deformation of the substrate reflects 
the Gaussian beam profile. The majority of the elastic 
strain energy is stored in the cavity substrate. Strong 
local deformations occur at the boundary between sub- 
strate and spacer extending into the spacer. At a distance 
exceeding this critical depth a homogeneous energy den- 
sity distribution is obtained. 



A. Local spacer deformations 

Following the discussion of Figure [2] it is clear that the 
analytic estimate for the spacer contribution according to 
Equation ^ cannot hold for the full spacer length as it 
neglects local deformations arising from the non-uniform 
pressure distribution on its front face. 

To illustrate the size of the effect, we estimate and sim- 
ulate the thermal noise following equations (6|8) for an 
all-ULE cavity with the dimensions and material param- 
eters given in Table [T] 
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FIG. 2: Deformation and contour plot of the elastic strain 
energy density in the mirror substrate and spacer for a Gaus- 
sian pressure profile with a 2 mm waist on the mirror surface. 
The color coding corresponds to the logarithm of the energy 
density pv in SI units. 



The results of the calculation are shown in Table HU 

The major part of the total strain energy caused by 
the static deformation is stored in the substrates. The 
contribution of the substrate is in good agreement with 
the analytic equations Q on the level of 1%. This re- 
sult reflects the good approximation of the mirror surface 
treated as an infinite half-space, i.e. the beam waist being 
much smaller than the mirror dimensions. The coating 
contribution was estimated by Equation Q for both the 
analytic calculation and the simulation. For the spacer 
geometry given in table [I] Equation ([8]) underestimates 
the energy stored in the spacer by approximately 30%. 
To shine light on this mismatch it is illustrative to cal- 
culate the energy distribution along the spacer length. 
The resulting graph for slices with a thickness of 1 mm 
is shown in Figure [3] 

The results have been normalized to the energy Uq 
contained in a slice of a homogeneously loaded spacer, as 
described by equation Q in section While the energy 
stored in the central part of the spacer matches the value 
predicted by equation Q, excess energy from additional 
deformations is stored close to the mirror. The pene- 
tration depth amounts to roughly 10 mm. The excess 
energy can be attributed to stress acting non-uniformly 
on the front face of the spacer as illustrated in Figure 
|2j The analytic estimate can be reproduced by applying 
a homogeneous force along the spacer's front face as ex- 
pected. The amount of excess energy on the spacer ends 
depends on the design parameters of the cavity. 

The dependence of the strain energy on the length of 
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coating 


analytic result 


1.04 


16.89 


0.126 


FEM calculation 


1.50 


16.99 


0.127 



TABLE I: Parameters used for the FEM simulation 



TABLE II: Comparison of strain energies in nJ for an all-ULE 
cavity 
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position along spacer axis (mm) 

FIG. 3: Simulated strain energy in slices of 1 mm thickness 
compared to a homogeneously loaded spacer (Co) (see Equa- 
tion |8|). The estimate by Numata et al. [T] (dashed line) as 
well as the analytic estimate according to Equation |8| (solid 
line) are shown as a reference. 
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FIG. 4: Strain energy as a function of spacer length for a 
spacer diameter of 32 mm. The estimate by Numata et al. 
PQ (dashed line) as well as the analytic estimate according to 
Equation Q (solid line) are shown as a reference. 
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FIG. 5: Elastic strain energy as a function of spacer diameter 
for a spacer length of 100 mm. The estimate by Numata et 
al. [I] (dashed line) as well as the analytic estimate according 
to Equation (JsJ) (solid line) are shown as a reference. 



tional to L the effect of the local deformations becomes 
negligible for long spacers. 

The dependence of the strain energy on the spacer di- 
ameter is illustrated in Figure [5] in comparison to the 
l/Rsp dependence of the analytic result (equation ^) 
and pp. All three curves show a decrease of the strain 
energy with increasing spacer diameter. The excess en- 
ergy increases from 0.45 nJ for a diameter of 32 mm as 
given in the example before and saturates at ^0.8 nJ for 
a spacer diameter approaching 100 mm. Note that for 
those diameters the spacer energy is totally dominated 
by the excess energy close to the mirrors ends. 

In summary, equation Q holds for the central region 
of the spacer while it fails close to the spacer ends. The 
impact of the spacer contribution on the total thermal 
noise budget of the cavity is small for all-ULE cavities, 
where the dominant part of the noise arises from the mir- 
rors. However, it can be strongly enhanced if the spacer's 
mechanical losses exceed the ones for the mirrors (see sec- 
tion 



III C ) or if nearly confocal cavities with large beam 



waists on the mirror surface lead to a reduction of the 
mirror thermal noise. In these cases a careful treatment 
of the spacer noise in a finite-element analysis is manda- 
tory. 



the cavity is shown in Figure [4] in comparison with the 
analytic estimates from equation ([8| and from pp. 

For cavities longer than the penetration depth of the 
local deformation of roughly 10 mm the simulation data 
are offset from the analytic estimate by a constant excess 
energy of approximately 0.45 nJ. Thus, once the elastic 
strain energy caused by the local deformations at the 
spacer's front faces has been determined for a given cav- 
ity geometry, the clastic energy of a cavity of a different 
length can be calculated according to equation ([8]). As 
the analytic estimate predicts the energy to be propor- 



B. Support structure 

To obtain fractional frequency instabilities below a y = 
10~ 15 , modern optical cavities employ vibration insensi- 
tive mountings [SJ 1141 1351 136| . In one widely used design 
for horizontal cavities, the resonators are supported at 
four positions close to the Airy-points by elastic materi- 
als pTOl HU [36l [37] . As the mechanical quality factor of 
rubber type material is relatively low the contribution of 
the support to the thermal noise might not be negligi- 
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Sketch of a cavity supported by four elastic support 



ble. To estimate the thermal noise contribution from an 
elastic support we investigate the geometry depicted in 
Figure [6] 

A spacer of length L is supported by four rubber sup- 
ports of thickness d sup , area A sup and shear modulus G sup 
at the Airy points positioned symmetrically at a distance 
L sup from the center. As the pads are located far from 
the spacer's front faces, local spacer deformations as dis- 
cussed in section [III A| can be neglected and we can treat 
the spacer as a homogeneously loaded bar. We again use 
the fluctuation-dissipation theorem in the form of equa- 
tion ([!]) to estimate the thermal noise contribution from 
the supports. We can estimate the energy stored in the 
support structure by assuming that the contact surfaces 
of support and spacer are displaced along the spacer axis 
direction according to the elastic expansion of the spacer. 
The support pads are placed at a distance L sup from the 
center of the spacer. Opposite forces F applied on the 
front faces of the spacer displace each rubber surface by 
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If we assume the support pads to be fixed to the support 
structure and have negligible action back on the deforma- 
tion of the spacer the resulting shear energy for a support 
pad is given by 
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area 


2 mm 
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TABLE III: Parameters for simulation of the cavity support 



tribution from four support points is then estimated by 
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The distance of the support from the center can be es- 
timated for an elastic bar by L sup = L/(2\/3) according 
to [35]. Consequently, the displacement noise S sup (f) 
for the support is proportional to L 2 resulting in a frac- 
tional displacement noise Sy P \f) — Sx (f)/L 2 inde- 
pendent of the cavity length. Therefore, care has to be 
taken especially for long spacers to avoid excessive influ- 
ence of the support on the total noise budget. 

For a vertical mounting of the cavity at the mid-plane 
(see e.g. (T3] [20]) the shear of the support pads acts 
transversal to the driving force deforming the spacer and 
the strain energy (equation (111) is reduced by a factor 
of (a sp • 2R C /L) 2 . 

For the horizontal cavity design with the parameters 
given in tables |III| and [I] we can estimate that the sup- 
port of 4 viton pads amounts to roughly 1% of the dis- 
placement noise Sx(/) = 8.95 • 10~ 35 // m 2 /Hz created 
in an ULE spacer. For spacer materials of higher me- 
chanical quality however the ratio (^sup/^sp can increase 
drastically and therefore in future cavity designs the im- 
portance of the mounting scheme will increase. For a 
cavity made of silicon for example Young's modulus in 
the [111] direction is 188 GPa [3S] and the mechanical 
loss angle of the material is <j> = 2 • 10 -8 [3D], about three 
orders of magnitude below the one for ULE (see Table |l|. 
Consequently, for a cavity made out of silicon the total 
thermal noise of the four support pads exceeds the noise 
of the spacer by a factor of ~ 4. Cavities of this type 
of material are currently limited entirely by the coating 
noise. However once low-loss coatings with sufficient re- 
flectivities [JT] are available this picture might change. 



Mixed material cavities 



In this approximation the shear energy is only affected 
by the shear strain and the displacement of the boundary 
is given by the spacer displacement itself. As a result a 
support with small area A sup , large thickness d sup and 
small shear modulus G sup will reduce the strain energy 
and therefore the thermal noise. The thermal noise con- 



In present cavities with ULE mirrors the major contri- 
bution to the cavity's thermal noise arises from their mir- 
ror substrates. Substrate materials of higher mechanical 
Q-factor like fused silica (FS) have been used [17l [19] to 
reduce this noise contribution. The thermal noise of FS- 
mirrors is already dominated by the noise of the coating 
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FIG. 7: Contributions to the thermal noise of the optical Spacer length L (mm) 

length S x (/) of a cavity with and without thermal expansion 

compensation ring. FIG. 8: Contributions to the thermal noise of the optical 

length S x (f) as function of the spacer length L. 



(compare Figure Un, However, fused silica shows a rela- 
tively large room temperature CTE of around 5 • 10 _ ' /K. 
The large CTE-difference in such a mixed material cav- 
ity leads to an unwanted lowering of the zero crossing 
temperature of the cavity's CTE in the order of a few 
10 K. To compensate for this effect, either special mirror 
configurations have been applied [35] or additional ULE 
rings have been optically contacted to the back surfaces 
of the FS mirrors [19] . 

We use the FEM model discussed above to check 
whether the additional ULE rings of the latter approach 
corrupt the low thermal noise of the FS mirrors. 

To the cavity model shown in Figure [T] additional 6 
mm thick ULE rings with a diameter of 25.4 mm and a 
central bore of 9 mm have been added. Figure [7] shows 
the FEM results of the thermal noise contributions of the 
different cavity components. We find that the additional 
thermal noise of the ULE ring is much smaller than the 
noise contribution of the mirror substrate or the mirror 
coating. The ULE rings enlarge the effective thickness 
of the mirrors which results in a smaller spacer deforma- 
tion and therefore in a smaller spacer excess energy. This 
reduction in the corresponding thermal noise is even big- 
ger than the additional thermal noise contribution of the 
rings itself and we have the paradoxical situation that 
a combined material cavity with additional ULE rings 
show less thermal noise than without these rings. 

While for all-ULE cavities the spacer contribution is 
negligible, in mixed material cavities the ULE spacer 
starts limiting the performance of the cavity. In gen- 
eral, following the discussion from section [TTJ, the noise 
contribution from the spacer is proportional to L and 
consequently a length L can be found where the spacer 
contribution equals the mirror contribution. This effect 
is illustrated for a fused silica type substrate in Figure 
[HJ For a cavity with the parameters given in Table [I] this 
length would correspond to 500 mm. 



IV. CONCLUSION 

Following Levin's "direct approach" we have analyzed 
the fractional length change of state-of-the-art optical 
cavities driven by Brownian motion thermal noise. While 
the thermal noise contribution created by the mirror is in 
good agreement with analytic results, we find discrepan- 
cies to the estimates given in pQ for the spacer contribu- 
tion. We give a revised analytic equation valid for long 
spacers and find an additional contribution because of 
additional deformations at the spacer ends close to the 
mirrors. This excess energy can exceed the energy de- 
posited by linear elastic deformation and can dominate 
the spacer contribution to the thermal noise for large 
spacer diameters and also for short spacer lengths. 

Another non-negligible source for thermal noise can 
arise from the support structure of the material. We 
have given estimate equations for a cavity mounted on 
soft support pads from below. 

Finally we have investigated the thermal noise of 
mixed-material cavities. We confirm that an additional 
compensator ring made of ULE (see [TH]) has no effect 
on the total cavity thermal noise. 

Current cavities are mostly limited by the coating 
noise. However this may change as low-loss coatings are 
being developed like microstructured mirrors [41] or coat- 
ings based on monocrystalline ALrGai-^As heterostruc- 
tures |42) . Additional reduction of thermal noise may by 
achieved by using large mode diameters or higher-order 
transversal modes [I3J . In these cases the results given in 
this publication enable further reduction of the thermal 
noise from the spacer, e.g. using thick mirrors, long cav- 
ities and optimized mounting. With these precautions 
instabilities below 10 -16 seem to be possible. 
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